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Periodic orbits are the central ingredients of modern semiclassical theories and corrections to these 
are generally non-classical in origin. We show here that for the class of generic polygonal billiards, 
the corrections are predominantly classical in origin owing to the contributions from closed almost- 
periodic (CAP) orbit families. Furthermore, CAP orbit families outnumber periodic families but 
have comparable weights. They are hence indispensable for semiclassical quantization. 



PACS number(s); 05.45.Mt, 05.45. Ac 

There exists an approximate dual relationship between 
the spectrum of quantum energy eigenvalues and the clas- 
sical length spectrum of periodic orbits and this forms 
the central theme of modern semiclassical theories. This 
duality was first discovered for the case of hyperbolic dy- 
namics where all periodic orbits are isolated and unsta- 
ble j|] and it was subsequently extended to the case of 
marginally stable systems where periodic orbits occur in 
families [gj. In particular, within the class of billiard 
systems (particle moving freely inside an enclosure and 
reflecting specularly from the walls), such a duality ex- 
ists for polygons which are marginally stable and where 
periodic orbits with even bounces occur in bands 

In general, there are other (weaker) non-classical con- 
tributions that make the relationship only approximate 
Q and must be included at finite energy. For special 
cases however (the tilted stadium billiard |12| and the 
truncated hyperbola billiard ||l^) there is a source of 
classical correction as well. The aim of this paper is to 
show that for an entire class of systems, corrections to 
the periodic orbit sum are predominantly classical in ori- 
gin and are due to closed almost-periodic orbits. Also, 
because they are more numerous and have weights com- 
parable to those of periodic orbit families, such orbits 
are indispensable at finite energies. First, however, we 
shall outline the key steps leading to the semiclassical 
trace formula where periodic orbits are the sole classical 
ingredients. 

A convenient starting point is the relation ^ 
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where G and Gg.c. refer respectively to the exact and 
semiclassical energy dependent propagator (Green's func- 
tion) and {En\ are the energy eigenvalues. The approx- 
imate propagator, Gg.c. is obtained from a fourier trans- 
form of the semiclassical time dependent propagator ||] 
and for a billiard, 
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where the sum runs over all orbits at energy E — k'^ 
between q and q' having length l{q,q') and /i is the as- 



sociated Maslov index. For convenience, we have chosen 
the mass m = 1/2 and h = 1. 

In the limit fc — > oo, the amplitude term in eq. (|^) 
varies slowly and can be regarded as a constant. The 
contribution of a particular orbit thus depends solely on 
the rapidity with which its action changes as q is varied. 
For periodic orbits, the action S{q, q) does not vary along 
the orbit. Further, if it occurs in a band, the action does 
not vary in the transverse direction either and the q- 



integration merely picks up the area, 
band. Thus 
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X cos{krlp - 7r/4) - ^ ^ cos{kr'lp>). (4) 

p' r' — l 

where pav is the average density of states and the sums 
over p and p' run over primitive families and (marginally 
stable) isolated orbits respectively having length Ip. 

For an isolated unstable periodic orbit on the other 
hand, the transverse direction leads to closed orbits with 
actions that vary depending on the stability of the peri- 
odic orbit and its contribution to the trace depends on 
the eigenvalues of the Jacobian matrix arising from a lin- 
earization of the transverse flow. In contrast, closed non- 
periodic orbits generally have negligible weight since their 
action varies rapidly with q. In case of the tilted stadium 
however, there exists a family of closed non-periodic 
orbits for which the variation of action across the family 
(bouncing between the straight edges) is small and its 
contribution can be of the same order as the bouncing- 
ball periodic orbit family in the zero-tilt stadium. Due to 
its close association with orbit families in straight-edged 
billiards, it is surprising to note that diffraction P-pJ is 
still considered the most significant source of correction 
in generic polygonal enclosures. While this is certainly 
true when the set of allowed momenta is small, generic 
polygons have additional classical contributions that are 
by far more important. 

To underscore this point, consider an arbitrary poly- 
gon Ti obtained by perturbing another arbitrary polygon 
T. The slight change in the shape of the enclosure results 
in a slight change in the quantal eigenenergies so that the 
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structure of the length spectrum, S{x) (the power spec- 
trum of p{k) = 2kp{E)), is largely preserved and there 
are only minor variations in peak heights (see fig. |^). 
However, the spectrum of periodic orbit lengths in T and 
Ti are radically different as we shall shortly demonstrate. 
There is thus an apparent paradox which cannot be re- 
solved by invoking diffraction since their contributions 
are 0{k^^) at best compared to the 0{k~^^^) con- 
tributions of geometric periodic families. 
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FIG. 1. Length spectrum S{x) of the equilateral and 
(100l7r/30 00, 9997r/3000) triangle (referred to as Tl) The 
perimeter in both cases is 1. The arrows mark the positions 
of orbits that are periodic in the equilateral triangle but are 
almost-periodic in Tl. The full and dashed lines correspond 
to the equilateral and Tl triangles respectively. In both cases, 
the first 1100 levels have been used to obtain S{x). 

The change in length spectrum of periodic orbits upon 
perturbation is best illustrated by comparing the equi- 
lateral and Tl triangles. As in case of all rational poly- 
gons, the invariant surface of Tl is two dimensional and 
topologically equivalent to a sphere with g holes where 
g = 1 + (A/'/2) ^^(m^ — l)/ni where {miTr/rii} are the 
internal angles of the triangle and Af is the least common 
multiple of {rii}. Thus for the Tl triangle, g ~ 1000 while 
for the equilateral triangle, g = 1- Note that the number 
of allowed momenta values is 2Af so that if Af is large, 
the probability that two segments of a trajectory have 
a small angle intersection is large. Thus, even though 
the boundary is only slightly perturbed, the structure of 
the invariant surface changes radically. It may thus be 
expected that the spectrum of periodic orbit lengths in 
the two systems is very different as well. In the integrable 
case, these invariant trajectories live on the torus and are 
labelled by the winding numbers (Afi,Af2) which count 
the number of times the orbit goes around the two ir- 
reducible circuits. In the non-integrable case, very little 
prior information is available and we shall analyze 
the situation to demonstrate that the symbol sequences 
of periodic orbits in the equilateral triangle do not nec- 
essarily lead to periodic orbits in Tl. 

For the triangle enclosures, we shall use the symbols 
{1,2,3} for the three sides [|6j. A trajectory can then 



be labelled by a string of symbols S1S2 . . ■ s„ where Si G 
{1, 2, 3}. Thus a sequence 1323 denotes a trajectory that 
reflects off sides 1, 3, 2 and 3 respectively. Let us denote 
by Ri,i = 1,3 the 2x2 reflection matrices of the three 
sides. These can be expressed in terms of the angle di 
between the outward normal ( hi ) to a side and the 
positive X-axis : 



/ - cos(26',) - sin(26',) 
' ' sin(26',) cos(26l,) 
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Thus, for the sequence 1323, the initial and final veloci- 
ties are related by 
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where the superscripts f{i) refer respectively to final (ini- 
tial) velocities v whose components are Vx and Vy. It is 
easy to verify that when the number of reflections is odd 
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where = 2(6ii -I- 6I3 • . . + 6'„) - 2(6(2 +9i + . 
while for even number of reflections [n even) 
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where <^e = 2(6*1 + 6*3 ... -h On-i) - 2(6*2 -f 6*4 -f . . . 6'„). 

Obviously, the initial and final velocities can be equal if 
the resultant reflection matrix RsiS2...s„ has a unit eigen- 
value. For even n (the case of bands or families), the 
eigenvalues are e^*'^" so that the condition for the exis- 
tence of a unit eigenvalue is 



mod(27r) 
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For odd n on the other hand, the product of the eigenval- 
ues A1A2 = 1. The eigenvector corresponding to a unit 
eigenvalue is (sin(</?o/2), — cos((^o/2)) so that if a real or- 
bit exists with the sequence siS2...s„, its initial and 
final velocities are equal. 

In the event that a sequence repeats itself (denoted by 
S1S2 ■ ■ ■ Sn) and there exists a unit eigenvalue of the resul- 
tant matrix RsiS2...s„, stability considerations guarantee 
that a periodic orbit exists . When n is odd, the orbit 
is isolated where as when n is even the orbit exists in an 
equi-action family. 

Not all sequences are however allowed. Further, not 
all repeating sequences guarantee the existence of peri- 
odic orbits due to eq. (||). For the Tl triangle, it is clear 
that the set of repeating sequence are the same as in the 
equilateral triangle for short orbits. Eq. (^) however does 
not allow all of them to be periodic. For instance, the se- 
quence 1323 results in a bouncing ball family of periodic 
orbits in the equilateral triangle. In the Tl triangle how- 
ever, the eigenvalues for this sequence are exp(±i7r/1500) 
so that there can be no periodic orbit with refiections 
from these sides. A sequence that is however allowed and 
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leads to periodic orbit families in both triangles is 123123 
( this is distinct from 123) since the periodicity condition 
(eq. ^ is automatically satisfied. In general then, for an 
arbitrary enclosure close to the equilateral triangle, an 
allowed sequence that repeats itself in the equilateral case 
can be a periodic family only when each symbol occurs as 
many times in even places as in the odd places. Thus, 
corresponding to the sequence 3231231231, there does 
not exist any periodic orbit in the Tl triangle while a 
periodic family exists in the equilateral case. 



We have thus verified that the periodic orbits in the Tl 
and equilateral triangles are indeed different even though 
short orbits follow the same sequence due to the proxim- 
ity of the two triangles. Note that this observation holds 
in general for any arbitrary enclosure T. Upon perturba- 
tion, orbits follow the same sequence but the periodicity 
condition will not be satisfied for sequences that are pe- 
riodic in T. According to eq. (||) therefore, the peak 
positions and heights in the length spectrum should dif- 
fer and we shall now show that the similarity in length 
spectrum observed in fig. |^ is due to contributions from 
closed almost-periodic orbit families in Tl. 




FIG. 2. The unfolded trajectory 3231 (marked by an ar- 
row) is produced by successive reflections of triangle I to 
produce copies II, III, IV and V. For the equilateral case, 
copy I and V have the same orientation and the trajectory 
is periodic. For Tl, the orientations differ slightly as shown 
schematically in the right. As a result the orbits are closed 
but non-periodic. 



Consider a symbol sequence that repeats itself and ex- 
ists in both the equilateral and the Tl triangles. Further, 
assume that corresponding to this sequence, there does 
not exist any periodic orbit in the Tl triangle while a 
periodic orbit family does exist in the equilateral case. 
Examples of these are the sequences 3231, 3231231231 
{Ip = 1.5275) and 2312312312312131 [Ip = 2.5166). In 
every such case, one can construct "unfolded" trajecto- 
ries (which are straight lines) by successive reflections of 
the triangle about the sides where the collision occurs. 
For instance (see fig. || ), unfolded trajectories for the se- 
quence 3231 can be created by first refiecting the triangle 
about side 3. The copy (II) so obtained is then reflected 
about side 2, the resultant copy (III) reflected about side 
3 and finally (copy IV) about side 1. For the equilateral 
triangle, the final copy (V) has the same orientation as 
the initial copy (I) so that any line joining corresponding 
points in the initial and final copies is an "unfolded peri- 
odic orbit" . In the Tl triangle however, the final copy dif- 
fers marginally in orientation from the initial copy so that 
any line joining corresponding points in the two can only 
be a closed almost-periodic orbit. Obviously, at every 
point q there exists such a closed orbit with this sequence 
so long as the line joining the corresponding points (in I 
and V) lies entirely within the copies generated through 
reflections. Two such orbits separated by q± are shown 
in fig. ^ (right). It is easy to see that the orbits dif- 
fer in length by an amount = q± tan(A0) ~ q± A9 
if Ad is small. Note that the above analysis holds for 
other almost-periodic closed orbits as well (such as the 
sequence 3231231231) and any arbitrary polygon, T. In 
each of these cases A^ ~ q± AO = q± (pe so that the 
length varies slowly if the orbit nearly closes in momen- 
tum. 



The correct trace formula for an arbitrary polygon T 
can be derived by noting that for a closed almost-periodic 
family, l{q±) = l{0) + q±(pe where l{0) = k is the length 
of the orbit in the centre of the band and q± varies from 
—Wi/2 to Wi/2 where Wi is the transverse extent of the 
band. Assuming that k is sufficiently large, the amplitude 
{l/l{q±)) can be treated as a constant (l/h) and the trace 
formula for finite k is then 

sm{kipi^^Wi/2) 
X cos{kli - TT/A) — —J-, — 

kife'Wi/2 

oo , 

-EEiJfc-^(^-'v)- (10) 

p' r' — l 

In eq. (0), the sum over i runs over closed almost- 
periodic and periodic orbit families and U is the (average) 
length of such a family. Note that as /c ^ oo, the con- 
tribution of almost-periodic orbits {ipi^^ ^ 0) vanishes as 
so that eq. ( p^ reduces to eq. (^). For de Broghe 
wavelength, A >> irwiipi^, however, the (ith) closed 
almost-periodic orbit family contributes with a weight 
comparable to that of periodic families (©(l/Zc^/^)) and 
hence assumes greater significance than diffraction [ pO[ . 
Interestingly, such orbits clearly show up in eigenfunc- 
tions |^l| and this has been referred to as "scarring by 
ghosts of periodic orbits" since such a periodic orbit ex- 
ists only in a neighbouring polygon. Thus a direct resolu- 
tion of the paradox lies in closed almost-periodic orbits. 

To emphasize the importance of the angle between the 
initial and final momentum (v^e), we compare the power 
spectrum of three different triangles, T1,T2 and T3 with 
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the equilateral triangle in figure^. For the sequence 3231, 
iff. is maximum for T3 and minimum for Tl so that peak 
heights at 0.57 and its repetitions should be closest to 
those of the equilateral triangle for Tl and farthest for 
T3. This can indeed be verified from fig. ^ 




FIG. 3. A comparison of the length spectrum for 
four different triangles : EQUI - equilateral, Tl - 
(1.00l7r/3,0.9997r/3,7r/3), T2 - (1.0l7r/3, 0.997r/3, 7r/3) and 
T3 - (1.015137r/3, 0.984877r/3, 7r/3). The arrows are at 0.577 
and 1.154 corresponding to the sequence 3231. In all cases, 
the first 1100 levels have been used to obtain S{x). Note that 
Tl is practically indistinguishable from the equilateral curve 
while T3 is farthest from EQUI. The corresponding values of 
ipe for the four cases are : EQUI - 0, Tl - 0.0006677r, T2 - 
0.00 6 6 677r and T3 - 0.0100877r. In contrast, the peak at a; = 1 
remains unchanged for all 4 triangles since it corresponds to 
a periodic orbit (123123). 

The contributions of CAP families diminish with en- 
ergy in accordance with eq. (|l0| ) and can be observed in 
the length spectrum. In order to distinguish this from 
the contribution of periodic families, we shall consider 
the power spectrum, G(x), of p(fc)/fe^/^ 



G{x) = 
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such that for a fixed kp — ka, the peak height of peri- 
odic families remains unaltered irrespective of fc/j. Fig. ^ 
show plots of G{x) for the T2 triangle using two different 
A:-intervals : (21,521) and (200,700). In both cases, the 
peak height remains unaltered at x — 1.0 correspond- 
ing to a periodic family. The peak at a; = 0.57 however 
diminishes in height as the interval shifts to a higher en- 
ergy. Also shown is a plot for the equilateral triangle 
which remains unchanged so long as kp — ka is fixed. 

Precise checks (without using any window function) 
between the observed and expected peak height at x — 
0.57 show that the value expected from eq. ( p^ ) is 11.3 
while the observed height is 9.6. Undoubtedly, there are 
other sources of corrections but the dominant contribu- 
tion at this value of x is due to the closed almost-periodic 
familiy. 



CD 



'EC!UI_21_521' 
T2_21_521' 
T2 200 700' 



0.6 



0.7 



0.8 „ 0.9 




FIG. 4. A comparison of G{x) for 2 energy ranges for the 
T2 triangle together with a typical plot for the equilateral 
triangle (marked EQUI_21_521 with = 21 and k/} = 521) 
when kjs — kct — 500. Note the diminishing peak heights for 
the range k £ (200,700). 

Finally, though the examples chosen are close to the 
(7r/3, 7r/3) triangle, we wish to reiterate that closed 
almost-periodic families contribute away from the neigh- 
bourhood of integrable enclosures as well. To see this, 
consider an arbitrary triangle T. In its immediate neigh- 
bourhood, there exists an infinity of triangles {T^'^}, 
each with a distinct periodic orbit spectrum but having 
the same symbol sequence for short trajectories. Assume 
now that there exists a periodic orbit corresponding to 
the sequence Sk for the triangle T^^^. Then, for all other 
triangles in its neighbourhood, this sequence contributes 
an amount (nearly) equal to the periodic orbit contribu- 
tion of T^^^ provided i^WiLpt^ « X. Thus corresponding 
to every periodic family in each of the triangles {T^'^}, 
there exists an almost-periodic family in the triangle T 
whose contribution is comparable to that of periodic or- 
bit families in these neighbouring triangles. 

To conclude, we have demonstrated that closed 
almost-periodic orbit families are more numerous and 
have weights comparable to that of periodic families in 
polygonal billiards. They are thus indispensable for the 
semiclassical quantization of generic polygons. 
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